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ABSTRACT 


This paper presents an analytical study on how to use 
the residual motor heat of a solid-propellant rocket for 
propulsion purposes in full vacuum, through the introduction 
of an inert quenching fluid. Three fluids were considered: 
water, propane, and Freon 12, with Delta Launch Vehicle 
third stage rocket motors (X-258 or FW-4). The possibility 
of using this quenching technique to provide fine velocity 
control was explored. Water was found to be superior to 
the other fluids under the conditions assumed, and the 
method appears to be feasible for velocity-trim within 6 
minutes after solid-propellant burnout. The study includes 
a computer program to predict the effect of introducing the 
quenching fluid at various times subsequent to the peak 
motor temperature. 
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PROPULSIVE POTENTIAL OF RESIDUAL ROCKET MOTOR HEAT 


by 

Joseph Ansell 

Goddard Space Flight Center 


INTRODUCTION 

Critical injection missions require accurate control of the total impulse generated by the stages 
of a launch vehicle, especially the final stage. If an error occurs in the terminal injection velocity, 
the velocity can be corrected by a fine-trim control system. If the final stage is a solid-propellant 
motor, one method of fine velocity control is to inject a cold, inert fluid into the hot, expended 
motor (see Figure 1). While the motor is cooling and losing its energy, the inert fluid gains heat 


energy from the motor. Allowing the fluid to 
vaporize and flow through the rocket nozzle 
provides an advantage by use of its residual 
propulsive potential. The velocity gain can be 
adjusted to ensure a fine control over final 
velocity. The fineness of the control is de- 
termined by choice of a velocity-sensing sys- 
tem and the system of valves that permit only 
the required amount of inert fluid to enter the 
hot motor shell. 

This approach also mitigates outgassing 
(Reference 1), which often contributes to an 
error in total impulse. Outgassing causes two 
problems after the solid main grain burnout. 
First, the total impulse of the solid-propellant 
stage may be greater or less than that pre- 
dicted on the basis of static tests (Reference 2). 
Second, any outgassing that occurs after sepa- 
ration of the motor from the payload may 
cause the final-stage motor to collide with the 
payload, possibly damaging it or changing its 
attitude or course. A method of quenching a 
hot motor with an inert fluid to eliminate out- 
gassing has been successfully used in ground 
testing (Reference 2). 



Figure 1— Schematic of a system for the injection of an 
inert fluid into a hot, solid-propellant motor at a con- 
stant injection pressure. 
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The purpose of this study is: (1) to find the weight penalty and velocity loss due to carrying 
an inert fluid velocity fine-trim system, (2) to find the potential velocity that can be recovered 
from using such a system, and (3) to compare the velocity loss with the velocity recovered. A 
trade-off comparison will show whether it is feasible to carry an inert-fluid, velocity-trim system 
on a space vehicle that requires close control of the injection velocity. 

Appendix A gives the calculations of comparative values for three selected substances; Ap- 
pendix B lists motor-temperature measurements; Appendix C lists required motor specific heats. 
Appendix D is a reproduction of a computer program to predict the effect of introducing the quench- 
ing fluid at various times. Appendix E gives the calculations of equilibrium conditions for the three 
substances. Appendix F is a symbol list. 


ANALYSIS 

A solid-fuel motor retains, after firing, an amount of heat energy Q 0 . Since no convection 
occurs in a vacuum, subtracting the amount of heat Q rad radiated from the motor’s surface gives 
the amount of heat energy available from the hot motor at the peak motor temperature. 

H = Q 0 ' Qrad (1) 

Q 0 - W 3 C p (T-T 0 ) . ( 2 ) 

The heat radiated from the motor's surface per unit area per unit time is 

Qrad = H T4 -Ti 4 ) (3) 


and the total heat radiated is: 


Qrad = «r(T«-T*)At . (4) 

Therefore, the heat remaining in the motor at time t after maximum motor temperature is reached 
is given by: 


H = W 3 C p (T-T 0 ) - «rA(T*-V)t 


(5) 


and, since Tj is very much smaller than T, 

H - W 3 C p ( T - T o) " €cr A T 4 t . 


( 6 ) 
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Differentiating with respect to time (in order to obtain the time decay of heat from the motor) gives 


dH/dt = - ecrAT 4 + (W 3 C p - 4€<xA T 3 t)dT/dt . 

The heat remaining in the motor is now available for heating the inert propellant fluid; this 
heat, prior to the fluid flow through the hot motor, is 

dH 

H = H + t 

H = [w 3 C p (T-T 0 )] - 2 [ecrA T 4 t ] + [(w 3 C p t -4ecrAT 3 t 2 )dT/dt] . 

The drop in temperature caused by the flow of inert fluid may be found by equating the heat 
lost by the motor to the heat gained by the inert fluid: 


(7) 

( 8 ) 
(9) 


W(At )C p6 (T B - T 0 ) + W(At )k + W(At )C p7 j T + dT/dt - t + -jjr (At ) - T„| - W 3 C p [- dt/dt - t - dT m /dt(At )] -(10) 


Solving for dt m /dt gives 


, _ W(At)C p6 (T B -T 0 )+W(At)\+W 3 C p dT/dt-t+W(At)C p7 T+W(At)C p7 dT/dt-t-W(At)C p7 T B 

dTm//dt " ~ -W(At) 2 C p? -W 3 C p (At) ' 

To determine the heat available in the motor at any time for a given flow of inert fluid into it, 


dT 

H = W 3 C p (T-T 0 ) - 2eoA T 4 t + (w 3 C p t - 4eaAT 3 t 2 )dT/dt + W 3 C p ( At ) • (12) 

For any given time of injection of fluid into the hot motor and at any flow rate, the time after 
maximum temperature when the available heat falls to zero may be determined from Equation 12. 
A computer program helped determine this value (see Appendix D). To find the amount of fluid 
that can be injected, multiply t(H = 0) by the fluid weight flow rate. 

m = t(H = 0)W . ( 13 ) 

Equation 13 assumes immediate (flash) vaporization of the quenching fluid on entry into the 
motor chamber. This mass of inert fluid will supply a total impulse 

= mI vac (14) 


and, since m and I decrease with time, I,„. is a maximum when t = 0. 

* vac ' cor 
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Velocity penalty and boost: 


w 1 + w 2 

AV 1 8 I vaco ln W, +W 3 

W 1 +W 2 +W 4 +W 5 

AV 2 = S I vaco ln Wj +W 3 +W 4 +W 5 
Velocity penalty = AV p = AVj - AV 2 

W 1 + w 3 +W 4 +W 5 

AV 3 (t) = velocity regained = gl vac (t) In ^ + Wj + w ^ + Wj (18) 

AV 3 is maximum at time t - 0 . 


(15) 

(16) 
(17) 


COMPARISON WITH THEORETICAL RESULTS 


NOZZLE CLOSURE 


RELIEF GROOVE 


EXIT 

CONE LINER 


An analysis of this type of quenching system was made with the X-258 and FW-4 Delta third- 
stage solid rocket motors (References 3, 4, and 5). Actually, given the appropriate data, the de- 
rived relationships can be applied to any solid-fuel motor to which a quenching fluid is to be ap- 
plied after burnout. For the X-258 (S/N RH-47) 
motor (Figure 2) using estimated values for 
specific heat, mass, maximum temperature, 
and motor surface area, the approximate values 
for incremental velocity increase have been 
calculated. In this study it is assumed that 
maximum motor temperature does not exceed 
chamber outer wall-temperature measurements. 
This gives conservative values, since the max- 
imum inner wall temperature is considerably 
higher than that estimated from outer wall- 
temperature measurements. 


IGNITER 

ASSEMBLY 


FWD 
ADAPTER 



FWD 

INSULATION 


FIBERGLASS 

CHAMBER 


AFT 
INSULATION 


ALUMINUM SKIRT 
ATTACHMENT RING 


Figure 2— Allegheny Ballistics Laboratory X-258 
B-l rocket motor. 


The following results were obtained by calculating propulsive potential for three different 
inert quenching fluids. Carbon dioxide has previously been used experimentally (Reference 2) in 
simulated altitude tests for quenching, and water has been proposed repeatedly. Theoretical 
studies have provided criteria for the selection of a cold-gas propulsion system (Reference 6). 
Water, propane, and Freon-12 were selected for this study; the choice was dictated by their mole 
cular weight, heat of vaporization, and specific impulse. 


Water as a Coolant and Propellant 

Water (References 7 and 8), with its low molecular weight and high specific impulse, is an 
excellent coolant and propellant in spite of its extremely high heat of vaporization. Also, it costs 
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nothing and is readily available (Reference 1). An analysis of a quenching system that uses water 
as the inert material with an X-258 (S/N RH-47) motor gave the following quantities for the derived 
relations (for constants and calculations, see Appendix A, part 1): 

H(t = 0) = 7800Btu. 

Figure 3 shows the heat radiation vs time from the motor. The maximum quenching fluid mass for 
water for the flow rates used was 

W 5 (t = 0) = 8.89 lbm for a flow rate of 0.6 lbm/sec. 

Figure 4 shows the effect of water flow for flow rates of 0.03 to 0.6 lbm/sec. The maximum amount 
of fluid permits the total impulse recovery of 



Figure 3— Heat radiation versus time. 


Figure 4— Available H 2 0 mass versus time. 
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from specific impulses calculated by the method of References 9, 10, and 11. Figure 5 shows the 
total impulse recoverable with time for various water-flow rates. The velocity penalty for an 
X-258 motor carrying water in addition to a 300-lbm payload is obtained from Equation 17. 


V, = 7362.80 ft/sec, 
v 2 = 7218.46 ft/sec, 


v = 144.34 ft/sec, 

P 

V 3 (t = 0) = 113.94 ft/sec, 

Figure 6 shows the recoverable velocity as 
a function of time and flow rate, with the great- 
est velocity recovered at the highest flow rate. 
The above values show that the maximum attain- 
able velocity occurs when no quenching system 
is carried. However, the recoverable-velocity 



TIME WATER INJECTED AFTER 
MEAN TEMPERATURE (sec) 

Figure 5— Total impulse available 
from H 2 0 versus time. 


no water carried 

water carried and not used, with a 
carrying structure equal to 20 percent 
of the fluid weight 

velocity penalty 

maximum recoverable velocity from 
the quenching fluid. 



Figure 6— Recoverable velocity versus 
time and flow rate. 
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gain when a quenching system is carried and used incurs a relatively small loss (31 ft/sec in 7300 
ft/sec or 0.4 percent of av); the control of total impulse error must be assessed against this loss 
when a mission requires very low injection-velocity error. 

Comparison of Water, Propane, and Freon-12 

Table 1 shows the values for the same calculations as water (above) for Freon-12 and propane. 
For constants and calculations, see Appendix A, parts 2 and 3. Of the three substances studied for 
this report, water proved to give the best velocity recovery; propane was a very close second. Both 
Freon-12 and propane had a higher specific impulse at the high initial temperature of the motor 
than is shown in Reference 8, a low temperature study. However, their specific -impulse values 
relative to each other agreed fairly well with the referenced literature. 

Table 1* 


Comparison of Three Different Materials for Theoretical Maximum Values at t = 0. 


Measured 

Quantity 

Water 

Propane 

Freon-12 

H(t = 0) (Btu) 

7800.00 

7800.00 

7800.00 

M(t = 0) (lbm) 

8.89 

17.72 

40.19 

I tot (t = 0) (lbf sec) 

1261.49 

2507.07 

2655.59 

AV 1 (ft/sec) 

7362.80 

7362.80 

7362.80 

AV 2 (ft/sec) 

7218.46 

7081.49 

6753.61 

AV p (ft/sec) 

144.34 

281.31 

609.19 

Av 3 (t = 0) (ft/sec) 

113.94 

222.62 

226.17 

AV 3 - V p (ft/sec) 

-30.40 

-58.69 

-363.02 


’The last line fAV, V 1 represents the net profit (+) or loss (~) in incremental velocity upon expending all the fluid. Zero corre- 
sponds to the break-even point. 


ANALYSIS OF SYSTEM 

In this study, it has been noted that the final attainable velocity depends mainly on the solid- 
propellant motor residual heat and the choice of a few important inert-fluid properties. When water 
is used as a quenching agent, its high heat of vaporization is a definite shortcoming. A substance 
with a low heat of vaporization gives more efficient use of the motor heat, but the molecular weight 
of the substance should not be high enough to give an unduly low I sp . Another factor required is 
high specific impulse for the available peak temperature. Appendix E lists theoretical rocket per- 
formances for equilibrium conditions, with these inert fluids. The specific heat of the motor shell 
must also be considered (see Appendix C). So must its mass; the higher the mass, the more energy 
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does it provide lor a given maximum temperature and specific heat, but the more velocity penalty 
does it incur. The at, and therefore the maximum temperature of the motor, should be as great as 
possible. Insulation such as the foil used on the X-258 (S/N RH-47) works well for two reasons: 
it allows a high maximum temperature, and it retains the heat longer than without insulation. Sur- 
face area and emissivity of the motor shell have a slight (almost negligible) effect on heat leakage 
after time t = 0; the smaller they are, the more slowly will the heat be radiated from the motor's 
surface. 

Finally, the time that elapses following maximum temperature is very important; if a quench- 
ing fluid is not used quickly enough, the system may not allow the recovery of any velocity. 


CONCLUSIONS AND RECOMMENDATIONS 

The results of this study show that although the velocity lost by carrying the weight of an inert 
quenching fluid cannot be fully regained, the penalty is fairly small for water (30 ft/sec), somewhat 
higher for propane (59 ft/sec) and considerably higher for Freon 12 (363 ft/sec). From a theoret- 
ical standpoint, the residual rocket motor heat constitutes a feasible source of energy for potential 
application to a propulsive system. For those launch conditions that require close control of the 
injection velocity, injecting an inert fluid into the hot post-burnout combustion chamber should be 
a way to adjust the final velocity until the motor has lost its heat from radiation. The results are 
considered conservative because, for a slight increase in either the chamber temperature or the 
heat capacity, there might be a velocity gain instead of a penalty. 

It would be advisable to follow this theoretical study by experimental measurements of: (1) the 
actual residual heat energy after firing, and (2) the effect of introducing an inert quenching fluid 
at various flow rates (using adequate thrust and pressure instrumentation). It should then be pos- 
sible to establish the efficiency of the propulsive system under realistic conditions, to determine 
the resolution in fine velocity control with existing hardware, and to generate the criteria for se- 
lecting competitive velocity-correcting techniques. 
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APPENDIX A 

CALCULATION OF COMPARATIVE VALUES FOR 
WATER, PROPANE, AND FREON-12, AS COOLANT AND PROPELLANT 


X-258 (S/N RH-47) motor: water as coolant and propellant 

(a) Constants: 

W 3 * 50 lbm 
C p * 0.30 Btu/lbm°R 
T % 600°F = 1060°R = 590°K 
T 0 * 80°F = 540 °R 
e % 0.7 

CT = 5.672 X 10~ 5 erg/°K 4 sec cm 2 
= 3.30 x 10’ 15 Btu/°R 4 sec in 2 
T, % 0°K 
A % 3700 in 2 

dT/dt * 0.25°R/sec (see Figure Al) 

C = 1 Btu/lbm°R 
PS 

C = 0.48 Btu/lbm°R 
\ = 540 Btu/lbm 

T b = 212°F - 672°R (this assumes that the injection pressure is constant at 1 
atmosphere) 

I vac = 142.1 lbf sec/lbm at t = 0 (t = 590°K, A e /A t = 53.2) 
dl vac = -0.0181 lbf/lbm (see Figure A2) 
w, % 300 lbm 
w 2 ~ 550 lbm 
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tankage, etc., 

When: 

W x = 0.03 lbm/sec then dT ra /dt = -1.47°R/sec 
w 2 = 0.1 lbm/sec then dT m /dt = -5.46°R/sec 
W 3 =0.3 lbm/sec then dT m /dt = -16.7°R/sec 
w 4 = 0.6 lbm/sec then dT m /dt = -33.4°R/sec 
I ^ 258 lbf sec/lbm 

vaco 

(b) Calculations: 

Total heat available: 

H(t = 0) = 50 lbm x 0.30 Btu/lbm°R x 520°R 
H (t = 0) = 7800 Btu (see other values in Figure 3) 

Total possible mass: 

m(t = 0, w = 0.6) = 8.89 lbm (see other values in Figure 4) 

Total impulse: 

I ( t = 0, w = 0.6) = 1261.49 lbf/ sec (see other values in Figure 5) 
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Velocity penalty and boost: 


! 

AVj = (32.174) (258) In 300 + 50 


AVj = 7362.80 ft/sec 


AV 


2 


= (32.174) (258) In 


300 +550 + 1.78 +8.89 
300+ 50+ 1.78 + 8.89 


AV 2 = 7218.46 ft/sec 
AV p = 144.34 ft/sec 


AVj (t = 0) = (32. 174) (142.1 -0.0181x 0.6) In 


300 +50 + 1.78 +8.89 
300 + 50 + 1.78 


AVj (t = 0) = 113.94 ft/sec 


2. X-258 (S/N RH-47) motor: Freon-12 as coolant and propellant 


(a) Constants: 

C p5 = 0.22 Btu/lbm°R 
C p6 =0.15 Btu/lbm°R 
\ = 60.5 Btu/lbm 

T b = 028 °C = -18.5°F = 441. 5°R (6.72 atmosphere at 80°F) 
I vac = 66.7 lbf sec/lbm at t = 0 (T = 590°K A = 53.2) 
dI vac = -0.0093 lbf/lbm 

When: 

V/ j = 0.03 lbm/sec then dT m /dt = -0.07°R/sec 
Wj = 0.1 lbm/sec then dT m /dt = -0.82°R/sec 
W 3 = 0.3 lbm/sec then dT m /dt = -2.94°R/sec 
W 4 = 0.6 lbm/sec then dT m /dt = -6.12°R/sec 

(b) Calculations: 


Total possible mass: 

rn(t = 0, W = 0.6) = 40.19 lbm 
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Total impulse: 

I tot (t = o, W = 0.6) = 2655.59 Msec 


Velocity penalty and boost: 


AV„ 


= (32.174) (258) In 


300 + 550 + (1.2 x 40. 19) 
300+ 50+ (1.2 x 40.19) 


AV 2 = 6753.61 ft/sec 


AV = 609.19 ft/sec 

P 


V 3 Ct 


= 0) = (32. 174) (66.7 - 0.0093 X 0.6) In 


300 + 50 + (1.2 x 40. 19) 
300 + 50 + (0.2x 40.19) 


V 3 (t = 0) = 226. 17 ft/sec 


3. X-258 (S/N RH-47) motor: Propane as coolant and propellant 

(a) Constants: 

C p5 = 0.58 Btu/lbm°R 
C p6 = 0.40 Btu/lbm°R 
.V = 146 Btu/lbm 

t b = -42°C = -43.5°F = 416. 5°R (9.72 atmosphere at 80°F) 
I vac = 142.2 lbf-sec/lbm at t = 0 (T = 590°K A e /a t = 53.2) 
dl vac /dt = -0.0239 lbf/lbm 
When: 

w t = 0.03 lbm/sec then dT m /dt = -0.57°R/sec 
w 2 = 0.1 lbm/sec then dT m /dt = -2.48°R/sec 
w 3 =0.3 lbm/sec then dT m /dt = -5.01°R/sec 
w 4 =0.6 lbm/sec then dT m /dt = -15.94°R/sec 

(b) Calculations: 
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Total usable mass: 

m(t = 0, w = 0.6) = 17.72 lbm 



Total impulse: 

l tot (t = 0, W = 0.6) = 2507.07 lbf/sec 

Velocity penalty and boost: 

300 + 550 + (1.2 x 17.72) 


AV, = (32.174) (258) In 


AV 2 = 7081.49 ft/sec 


V = 281.31 ft/sec 

P 


300+ 50+ (1.2 x 17.72) 


v 3 (t = 0) = (32. 174) (142. 2 - 0.0239 X 0.6) In 
V 3 (t = 0) = 222.62 ft/sec 


300 + 50 + (1.2 x 17.72) 
300 + 50 + (0.2 x 17.72) 
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APPENDIX B 


MEASUREMENT OF MAXIMUM MOTOR TEMPERATURE AFTER FIRING 


Maximum temperatures were measured at nine separate points around the motor shell during 
and after firing (Reference 3). The temperatures at all points were averaged at different times 
after ignition to give the time of maximum motor temperature and the temperature average at that 
time. (See Table Bl.) 


Table Bl 

Averages of Maximum Motor 
Temperatures* of the X-258 RH-47 
Motor at Different Times After 
Ignition (Reference 3). 


t, sec** 

T, °F 

200 

576 

250 

596 

300 

603 

350 

602 

400 

599 


* Maximum average motor temperature: Approxi- 
mately 600°F at 300 seconds after ignition. 
Initial motor temperature - 80°F 
T = 600°F - 80°F = 520°F 
**Note: This time should not be confused with the 
t used in developing the thermodynamic equa- 
tions. That t is the number of seconds past max- 
imum motor temperature. This t - 300 + x. 


Change in temperature (AT) is the difference between 
initial motor temperature and maximum temperature. 

Table B2 gives a comparison of maximum temperatures and 
AT ? s for different motors and insulation. Figure A1 gives 
values for the time decrease of heat from the motor's sur- 
face, plotted against temperature starting with a fixed tem- 
perature of 360°F. 


Table B2 


Maximum Motor Temperatures of Different Type Motors. 



Motor Type 

Measured 

Quantity 

X-258 

(References 3 and 4) 

FW-4 (Reference 5) 

Insulation 

Foil 

No Foil 

No Foil 

Max T, °F 

570°, 600° 

425°, 440° 

340°, 390° 

AT, °F 

490°, 520° 

345°, 360° 

265° , 315° 
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APPENDIX C 


EFFECT OF SPECIFIC HEAT OF MOTOR SHELL ON AV 3 FOR WATER 


As can be seen from Appendix A, part 1, 
Figure 6, and Table Cl, the specific heat C p of 
0.30 which was estimated for the motor shell 



0.2 0.3 0.4 0.5 

C p ( Bhj/ Ibm °R ) 


Figure Cl— AV versus C D at t = 0. 


Table Cl 


C p T s of Motor Shell and Masses at H 2 0 
Required for Certain Velocities. 


V 3# ft/sec 

H 2 0 Mass Required 
(lbm) 

Necessary Motor C 
(Btu/lbm °R) P 

113.94 

8.89 

0.30 

144.34 

11.35 

0.37 

160.00 

12.33 

0.41 

180.00 

14.11 

0.46 

200.00 

15.89 

0.52 


will not provide the motor with enough heat 
energy to allow enough water to be used to re- 
coup the velocity penalty for carrying the water, 
A separate study gives the relation between in- 
crease in c p , increase in fluid mass, and con- 
sequently increase in Av 3 . The data from this 
study can be seen in Table Cl, and Figure Cl. 

It shows that to return just the velocity penalty 
at time t = 0, a motor shell with a c of at 
least 0.37 Btn/Ibm °R must be used. 
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APPENDIX D 
PROGRAM LISTING 


C WATER «CCKET 

C ENERGY REMAINING IN HOT MOTOR 

C TOTAL MASS, IMPULSE? AND INCREMENTAL VELOCITY AVAILABLE 

C HIGH TEMPERATURE PROGRAM 

C FOR THOSE MATERIALS WHOSE BOILING POINTS APE ABCVE STARTING TEMPEPATURF 
D THEN SIGN F < 14 ) , AM. ( 14 ) , T I MP ( 1 4 ) , VE L ( 14 ) , T ( 1 4 ) 

C I MENS ION DEL T( 1A, I A ) , TT ( 1 A) , T TD( 14) ,HH (1A ) ,HDUM ( 1A ,1A ) 

«FAP< 5,76) NCASL 

76 FORMAT (IE) 

CO ? 1=1, NC A SE 

PFAD( =>77) SIG»CM,CP*TEM,TEMC,EP»ARrA,DTCT 

77 FORMAT ( 8( F1C .E ) ) 

READ! 3, 7.° ) D I OT , DT WDT , TEMIC »C P 1 ,CP2 , AL AM , V AC l 0 ,G 
7 5 FORMAT ( 8 ( F 1 C . A ) ) 

READ! 5 » b C ) FLOW 
80 FORMAT! F 1C. 3 ) 

WRITE! 6 * 75) L 

79 FCRMAT113HCASC JUMLGR * I 5/7 / 1 1 CH TIME EELTA 

1 T ENERGY MASS IMPULSE VELECITY 

: TEST H////) 

EG 1 1=1,13 
T ( 1 )=-2C . 

CO 2 J = i, 13 
T ( J ♦ I ) =T ( J >+2C. 

CC=I 

CELT! I, J ) =T ( J + l )-( ( 3C .*01) >-30. ) 

IF ( DFLT< I, J ) J 3, 1C, 1C 
A EEL T I I,J ) =0 . 

10 F 1 = ( CM *CP*( TEM- TEMQ ) ) 

F2 = l 2.*EPAS IC*AREA*( TEM**4)*T I J+i > ) 

F 3 = 1 EM*CP*T( J+l X--OTOT) 

FA = (4.«EP*Sl G*AREA* ( TEMA*3) ADTDTA ( T ( J + l ) ** 2 ) ) 

F5 = ( OM*CP*DT»,OT*DELT ( I , J) ) 

F C U M ( I,J )=H1-F2+H3-H4+H5 
IF { F*CUM ( I, J ) ) 13,13,2 
15 v V = 32 . 

TT( I)=T< J + l ) 

TTIH I ) = C E L I ( l , j i 
18 VV = V V / 2 . 

F 1 = ( CM *CP A ( TEH- T EM 0 ) ) 

F2 = < 2 .*EPAS IGAAREA* ( TEMA* 4 ) $ TT ( I ) 

F 3 = ( E M A C P * T T ( I ) * 0 T 0 T ) 

F4*( 4 ,*rp*s IG*AREA* ( TEM**3)*0TDT* ( TT I I ) **2) ) 

F 5 = ( FM ACPADTWDTATTDI I ) ) 

FF( I )=F 1-F2 + H3-H4 + H5 
IF ( HF ( I ) ) 30,17,32 
30 IT ( HF ( I l + .Cl) 31,17,17 
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21 T T { I ) = T T ( I >-VV 
TTD( I ) = T T C ( I ) — V V 

GO TO 13 

32 IF ( HE { I )- • C 1 ) 17, 17,33 

33 TT ( I ) = T T { I )+VV 
TTDI I ) = T T D ( I )+VV 

CC TO 12 

17 AMI [ ) = F L G W * T T 0 { I ) 

8 H I ) = AN* ( I )*{CP1*< TEMK.-TEMO+ALAM+CP2* (TE M+( CTCT*TT I I ) ) + ( CTWCT*TTC t 
1 I > ) - T E M B > ) 

9 TIMPl I )=AM( I )*( VAC IC+( D 10 T* TT I I ) ) ) 

V1 = ALCG( < 35C. + I 1.2*AM I ) ) )/<35C. + ( .2*AM ! ) > >) 

VEL(I)=G*( VACI0 + (DIDT*TT(I) ) >*V1 
CO TO 1 

2 CONTINUE 
1 CONTINUE 

WR I TE ( £, 8 1 ) I TT( I ) ,TTD(I) ,Hl I ) tAK(l) ,TIMP U) ,VEL ( I) »FH( I ) , 1=1, 13) 

81 FORMAT ( 5X,7F15.3/) 

3 CONTINUE 
55 STOP 

END 


TER M DEFINITIONS 

F ENERGY REMAINING IN THE HCT MOT C R 

AM MASS OF WATER THAT CAN BE ADDFD TC THE HCT MGTCR 

T I M P IMPULSE AVAILABLE FROM THE WATER 

V EL ACCITIQNAL VELOCITY AVAILABLE FRCM TEE WATER 

T TIME ELAPSED SINCE MAXIMUM MOTOR TEMPERATURE 

CELT TIME ELAPSED FOR WATER FLCW 

SIG ST E PEAN-BOL T 2 MANN CONSTANT 

EM TCTAl MOTOR MASS 

CP AVERAGE MOTOR SPECIFIC HEAT 

TEM MAXIMUM MOTOR TEMPERATURE 

TEMO INITIAL MOTOR TEMPERATURE 

EP MCTCR SHELL EMISSIVITY 

AREA MCTCR SURFACE AREA 

CTCT MCTCR TEMPERATURE LOSS PER TIME TC SFACE 

DIET IMPULSE LOSS PER TEMPERATURE FRCM TEE WATER 

CTWDT MOTOR TEMPERATURE LOSS PER TIME TC WATER 

TEMO FOILING POINT OF WATER AT ATMOSPHERIC PRESSURE 

CPI SPECIFIC HEAT OF WATER LIQUID 

CP 2 SPECIFIC HEAT OF WATER VAPOR 

A L AM FEAT OF VAPORIZATION FOR WATER 

VACIC VACUUM SPECIFIC IMPULSE OF WATER AT THF MAXIMUM MOTCP TEMPERATURE 

G GRAVITY CONSTANT 
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APPENDIX E 


EQUILIBRIUM CALCULATIONS FOR 
WATER, PROPANE, AND FREON-12 



THEORE' ICAl rocket performance assuming equilibrium composition during expansion 

FROM AN ASSIGNED TEMPERATURE AND AT VARIOUS CHAMBER PRESSURES 
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ADDITIONAL PRODUCTS WHICH WERE CONSIDERED BUT WHOSE MOLE FRACTIONS WERE LESS THAN 0.000005 FOR ALL ASSIGNED CONDITIONS 



THEORETICAL ROCKET PERFORMANCE ASSUMING EQUILIBRIUM COMPOSITION DURING EXPANSION 
FROM AN ASSIGNED TEMPERATURE AND AT VARIOUS CHAMBER PRESSURES 
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ADDITIONAL PRODUCTS WHICH WERE CONSIDERED BUT WHOSE MOLE FRACTIONS WERE LESS THAN 0,000005 FOR ALL ASSIGNED CONDITIONS 


THEORETICAL ROCKET PERFORMANCE ASSUMING EQUILIBRIUM COMPOSITION DURING EXPANSION 
FROM AN ASSIGNED TEMPERATURE AND AT VARIOUS CHAMBER PRESSURES 
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ADDITIONAL PRODUCTS WHICH WERE CONSIDERED BUT WHOSE MOLE FRACTIONS WERE LESS THAN 0,300005 FOR ALL ASSIGNED CONDITIONS 
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ADDITIONAL PRODUCTS WHICH HERE CONSIDERED BUT HHOSE MOLE FRACTIONS HERE LESS THAN 0.000005 FOR ALL ASSIGNED CONDITIONS 
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APPENDIX F 
SYMBOL LIST 

A = surface area of solid-propellant rocket motor 

C p = average specific heat of motor = (M x C pl + m 2 C p2 +m 3 C p3 + M 4 C p4 + M 5 C pS yw 3 
c = specific heat of outer wall of motor 

C p2 = specific heat of metal parts (ring and bolts attached to motor) 

C p3 = specific heat of graphite throat 
C p4 = specific heat of motor liner 
C p5 = specific heat of exit cone 
C p6 = specific heat of inert fluid 
C p7 = specific heat of inert fluid vapor 
g = gravitational constant 

H = total heat in motor at maximum temperature 
H = available heat from motor at any time 
l tot - total impulse available from inert fluid 

I vac = vacuum specific impulse of inert fluid at the maximum motor temperature 
I vaco = vacuum specific impulse for solid-propellant rocket motor at t = 0 
M, = weight of outer wall of motor 

M 2 = weight of metal parts (ring and bolts attached to motor) 

M 3 = weight of graphite throat 
M 4 = weight of inner motor liner 
M 5 = weight of exit cone 

m t = weight of fluid that can be injected into hot motor at any time 
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m 0 = weight of fluid actually injected into hot motor 

Q 0 = energy generated and retained by rocket motor after main grain burnout 
Q rad = energy radiated by motor’s surface 
°R = degrees Rankine 

T = maximum temperature of motor shell after firing 
T 0 = temperature of motor shell before firing 

T t = temperature of space outside motor 

T b = boiling point of inert fluid at injection pressure 

t = time elapsed after maximum temperature (T) is reached 
At = time elapsed for flow of fluid into hot motor 

AV t = velocity increment provided by solid-propellant motor without quenching system 

av 2 = velocity increment provided by solid-propellant motor when carrying the injection sys- 
tem but not using it 

AV 3 = incremental velocity increase when using the system at any time 
AV p = V, - V 2 = velocity penalty 
W = weight flow rate of inert fluid into hot motor 
Wj = weight of payload 
w 2 = weight of motor 

w 3 = total weight of empty solid propellant motor = M, +m 2 +m 3 +M,, +M S 
w 4 = weight of inert-fluid-carrying system without fluid 
w s = weight of inert fluid to be used 
dH/dt = loss of heat from motor with time 

dl/dt = loss of specific impulse of inert fluid with time as temperature drops 
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dT/dT 

d T m / dt 

€ 

cr 

\ 


temperature decay rate from hot motor to space surrounding motor 

temperature decay rate when inert fluid flows 

emissivity of motor shell 

Stephan-Boltzmann constant 

heat of vaporization of inert fluid 



